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NUMERICAL ANALYSIS OF CONTINUOUS 
FRAMES IN SPACE 


By JAMES MICHALOS,' M. ASCE 


1 Prof. of Civ. Eng., Iowa State College, Ames, Iowa. 


SYNOPSIS 


A numerical procedure, involving successive corrections, is presented for 
determining moments in continuous girders and frames whose members join at 
arbitrary angles in space. The component parts of these structures may be of 
constant or variable cross section. 


INTRODUCTION 


An analytical investigation of the structural behavior of a continuous build- 
ing frame is usually made by neglecting the torsional resistance of girders at 
right angles to the particular bent under study. Usually, this simplification 
has not caused great concern among designers. The effect of these girders, and 
the effect on these girders, can be appraised, providing one knows—or can ap- 
proximate—the relationship between angle of twist and torque for such mem- 
bers. An analysis that includes the torsional resistance of members at right 
angles to a frame can be made by application of the ordinary method of moment 
distribution? to all members framing into each joint. This has been explained 


2‘ Analysis of Continuous Frames by Distributing Fixed-End Moments,” by Hardy Cross, Transactions, 
ASCE, Vol. 97, 1932, 


3*'Design of Reinforced Concrete in Torsion,” by Paul Andersen, [bi/., Vol 103, 1938, 1503 

4" Reinforced Concrete Structures,” by Dean Peabody, Jr., Jolin Wiley & Sons, Inc., New York, N. Y., 
2d Ed., 1946, p. 419. ‘ 

Analysis of Three-Dimensionr,) Beam and Girder Framing," by "hil M. Ferguson, Journal, A.CT., 
Vol. 22, September, 1950, p. 61. 


Structures may be built whose members frame into each other at various 
angles in space rather than at right angles. It seems probable that require- 
ments of clearance, space, and architecture wil! accelerate the use of structures 
having nonorthogonal framing, not only in building construction but also in 
connection with urban expressway construction. Simplifications are neither as 
obvious nor as acceptable as in those cases in which only relatively small 


261-1 


+ 
= 
4; 


moments are introduced into the joint by members at right angles to a bent. 

In this paper a numerical procedure is presented for analyzing girders and 
frames that are continuous in space. It is limited herein to straight members, 
but these may be of constant or variable section. Furthermore, the members 
may be oriented in space in any arbitrary fashion, both with respect to their 
longitudinal axes and the principal axes of their cross sections. In general, the 
paper deals with structures in which there is no joint translation, or in which 
such movement is negligible, but the handling of cases in which translation 
occurs is outlined and illustrated by means of a numerical example. 

The method involves successive corrections and can be carried to any de- 
sired precision. It is an extension into space of the philosophy of moment dis- 
tribution, and it is based on a concept wherein all defined properties and mo- 
ments are referred to a system of rectangular axes. 

Expressions for the physical constants used in the numerical procedure 
are derived for the general case of a member in space subjected to mo- 
ments about each of three rectangular axes, none of which is parallel to the 
longitudinal axis of the member or to the principal axes of the cross section. 
Simplified expressions are obtained for several special, or more usual, cases. 
The numerical procedure is then outlined, and examples are presented. 


STIFFNESS AND CarRRy-Over Factors 


When the method of moment distribution is used in the analysis of a frame 
lying in a plane, the unbalanced moment at each joint is distributed to all mem- 
bers framing into that joint. The portion of the unbalanced moment taken by 
each member is in proportion to its “‘stiffness,’”’ which is defined as the moment 
necessary to rotate one end of a member through a unit angle when the other 
end is fixed‘and there is no relative displacement of the ends of the member. 
But a distributed moment, resulting from rotation at one end of a member, 
induces moment at the other (fixed) end. This moment is determined by 
multiplying the distributed moment by a “carry-over factor,”’ which is defined 
as the ratio of the moment induced at a fixed end to the moment producing ro- 
tation at the other end. , 

In the method presented in this paper, all moments are referred to a system 
of three rectangular axes, and a joint is considered free to rotate about one axis 
atatime. Consider a member that is not parallel to any of the three orthogonal 
planes associated with an arbitrarily chosen system of rectangular axes. If 
moment is applied about one of the reference axes at one end—and if rotation is 
prevented about the other two axes at that end and all three reference axes at 
the other end—carry-over moments will be induced at each end about the 
axes of no retation. Stiffness factors will be required for rotation about each of 
the three reference axes. Five carry-over factors will be required for each of the 
three axes; two for the end at which rotation has occurred, and three for the 
opposite end. The adopted sign convention is shown in Fig. 1, wherein positive 
end moments and rotations are represented vectorially. The reader should use 
the right-hand rule for the sense of moments and rotations. Thus, if the thumb 
of the right hand is pointed in the direction of the arrow, the curved fingers 
indicate the direction of the moment or rotation. 
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rotations 


MOMENTS (b) ROTATIONS 


Fie. 1.—Postrtve Moments aNnD ROTATIONS AT THE ENDS OF A MEMBER 


Several views of a member supported at both ends are shown in Fig. 2. In 
addition to an arbitrarily chosen set of rectangular axes, z, y, and z, an auxiliary 
pair of axes, z, and 2, is introduced, but only as an aid in the derivation of 


(co) PROJECTION ON xz PLANE (0) PROJECTION ON xz PLANE 


x, 
8 


y 


(b) TRUE LONGITUDINAL VIEW (b) TRUE LONGITUDINAL VIEW 


(c) TRUE ENO VIEW AT A (6) SECTION AT B {c) TRUE END VIEW AT A 


Fie. 2.—Moment Fig. 3.—Roration 
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expressions for stiffness and carry-over. They have no role in the procedure 
itself. The z)-axis is parallel to the rz-plane, and the z,-axis forms a right angle 

vith the longitudinal axis of the member, its projection in the zz-plane coin- 
ciding with the projection of the longitudinal axis of the member in the same 
plane. The axis 1-1 and the axis 2-2 are the principal axes of the cross section. 
The angles y and ¢ are true angles, and a is the projection of an angle upon the 
xz-plane. The positive sense of these angles is indicated by arrow. 

If a moment K, applied about the z-axis at end A of the member in Fig. 
2(a) produces a unit rotation about the same axis, while no rotation is allowed 
to take place about the z-axis and the y-axis at end A or about any of the axes 
at end B, the moments necessary to prevent rotation about these axes can be 
written as shown in Fig. 2(a). The moment K, is the defined stiffness about 
the z-axis at end A, and C,, and C,, are carry-over factors to the z-axis and the 
y-axis, respectively, at the same end, whereas C’,,, C’,,, and C’,, are carry-over 
factors to the far end, B. The contributions of the moments K,, C,,, K,, and 
C,. K, to torsional moment and to bending moments about principal axes at 
end A, as shown in Figs. 2(b) and 2(c), can be obtained as follows: 


M, = K, sinasint — C,, K,cosasin{ + C,, K, cost 

Mn = K,cosacosy — K,sinacosfsiny + C,, K,sin a cosy 
+ C,,K,cosacos{siny + C,, K,sin¢g siny 

Mi. = K, cosasiny + K,sinacos{ cosy + C,, K,sinasiny 
— C,,K,cosacos{ cosy — C,, K,sin{ cosy 


The bending moments induced about the principal axes at B are C'"4¢ Mai 
and C*-24— Mas, as shown in Fig. 2(d). The terms C'~'4z and C?-*,4p are 
carry-over factors, from A to B, with respect to axis 1-1 and axis 2-2, respec- 
tively. The value of these carry-over factors is 0.5 for prismatic members. 
Values for nonprismatic members can be obtained from available charts or 
tables.® 

* “Handbook of Frame Constants," Portland Cement Assn., Chicago, Ill., 1947. 

In Fig. 3(a) the end A of the member has been rotated through a unit angle 
about the z-axis only. The contributions of the unit angle to twist and to 
rotation about principal axes at end A, as shown in Figs. 3(b) and 3(c), can be 
obtained as follows: 

6, = sinasing 
6.1 = cosacosy — sinacos¢ siny 
6.2 = cosasiny + sin acosf cosy 


This angle of twist and these rotation angles can be related to the moments 
of Eq. 1 by means of stiffnesses kai, kn2, and The terms k,; and ky,» represent 
the bending stiffnesses associated with moment distribution. Their values, for 


4 ET, 


prismatic members, are L and 7 , respectively. In these expressions, 
4 


E is the modulus of elasticity in tension and compression, L is the length of the 


member, and J, and J, are the moments of inertia about axis l—] and axis 2-2, 
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respectively. For nonprismatic members, values of ka: and kn: are available 
elsewhere.* The term k, represents the moment necessary to twist one end of a 
member through a unit angle when the other end is fixed. Its value for pris- 


matic members is Gd G being the modulus of elasticity in shear and J being a 


L ’ 
torsion factor, such that G J represents the torque necessary to twist a member 
through a unit angle per unit length. The value of k, for nonprismatic members 
can be determined through the use of the shear and torsion analogy.’ 


7Laterally Loaded Plane Structures and Structures Curved in Space,"’ by Frank Baron and James 
P. Michalos, Transactions, ASCE, Vol. 117, 1952, p. 279. 


Three equations, relating end moments and end rotations about the longi- 
tudinal axis and the principal axes, can be written as follows: 


(sina sing K, — cosasingé C,,K,+ cos{C,, K.) = sinasin¢.. (3a) 

t 

(cos a cosy — sinacos{siny) K, + (sina cosy + cos a cost 
ni 


X siny) C., K, + (sing siny) Cy, K,] = cosacosy — sina cos¢ siny. . (3b) 


[(cos asiny + sin acos{ cosy) K, + (sinasiny 
n2 
— cos a cos { cosy) C,, K, — (sing cosy) C,, K,] 
= cosasiny + sina cosf cosy. . (3c) 


In an identical manner, similar groups of equations are obtained as the 
result of applying a moment K, about the z-axis and a moment K, about the 
y-axis in Fig. 2(a). 

From the application of K,— 


(cos asing K, — sinasing C,, K, — cos{C,, Kz) = cosasin¢..{4a) 
t 


[(sin a cosy + cos asiny cos ¢) K, + (cosacosy 
ni 


— sinasiny cosf{) K, + sinpsin¢ C,, Kz] 
= sinacosy + cosasiny cos¢.. (4b) 
and 


{(sinasiny — cosacos ycos ft) K, + (cosasiny 


+ sin a cosy cos C,. K, — cospsin{¢ C,, K,] 
= sinasiny — cosacosy . (4c) 
From the application of Ky— 


¢ K, + sin asing C,, Ky — cosasin{t C,, Ky) = cosf.... (5a) 


[siny sint K, + (cos acosy — sinasiny cost) C., Ky, 


+ (sin a cosy + cos asiny cos ¢) C,, K,] = sinysin¢. . (5b) 
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a 
and 


[— K, + (cosasiny + sina cosy cos C,, K, 
+ (sinasiny — cos a cosy cos C,, K,] =— (5c) 


From the solution of the three groups of simultaneous equations (Eqs. 3, 4, 
and 5) expressions for the three stiffnesses (K,, K,, and K,) and the six near-end 
carry-over factors (C,,, Cy., Cys, Ca, Cy, and C,,) were obtained and are shown 
at the top of Fig.4. The coefficients, NV, have been introduced for convenience, 
and they are defined in Fig. 5. These coefficients are functions of the angles 
only, and some of the expressions are greatly simplified—as shown in the bottom 
half of Fig. 5—in the usual case in which the angle y is equal to zero. 

The carry-over moments at end B of Fig. 2(a) can be written in terms of the 
torsional and bending moments at the same end as follows: 


=— M,sinasing + C'"42M,, (cosy cos a — siny cost 
x sin a) + (siny cos a + cosy cos ¢ sin a) 
C’",.M, = — M, cost + sing — sin¢ >... . (6) 


C’,,.M, = cos a + C'"4gM,, (cosy sin a + siny cost 
cosa + C2 sin a — cosy cos cos @) 

; 1e right-hand sides of Eqs. 1 are substituted for M,, M,:, and M,» in 
ae #3 Eqs. 6, expressions for the far-end carry-over factors can be obtained in terms 

7 a a of the two carry-over factors at A and the two carry-over factors that are re- 
AS ad: lated to the principal axes. The resulting expressions for C’,,, C’.,, and C’,, are 

°. given in Fig. 4, in terms of coefficients, N, which are defined in Fig. 5. 


z 

+ Ky = ke + — Kneis Nis) Cyz Cay = Cyz 

Ky = ky + Kn, + Cyn * Cuy* 

y 


Ch, Ny ~ + (Ne + + Nig) + Gag (Ny +Cyz Nig) 

Cyn + + ) + + Cya Neo) 

-CayNg -CayNig) + CryNz ) + +CryNo ) 


Fig. 4.—General Expressions for Stiffness and Carry-Over Factors 
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Similarly, the application of a moment K,, instead of K,, at A in Fig. 2(a) 
would lead to expressions for the three far-end carry-over factors C’,,, C’,2, and 
yz, and the application of a moment K, at A would lead to expressions for C’,,, 
C’,,, and C’,,. The expressions for these carry-over factors are included in Fig. 
5. In each case, the first subscript indicates the axis to which moment is 
carried, and the second subscript indicates the axis about which moment is 
applied. 

The expressions for stiffness and for carry-over factors, as shown in Fig. 4, 
are for the general case. However, for the majority of cases which may be en- 
countered in practice, the expressions are much simplified, as shown in the 
bottom half of Fig. 5 and in Fig. 6. 


THe NUMERICAL PROCEDURE 


Before the numerical procedure can be applied, it is necessary to compute 
the stiffnesses and the carry-over factors. In the great majority of cases, y is 
equal to zero, and the expressions are greatly simplified. If, in addition, ¢ is 
equal to 90° or to zero, the expressions in Fig. 6 are directly applicable. 

Only relative values of kn, kn2, and k; are necessary, and for prismatic mem- 
bers both C'~'4, and C?-24g are equal to 0.5, thus further simplifying the ex- 
pressions for carry-over factors. In dealing with prismatic members, or with 
symmetrically haunched members, the required stiffness values are the same for 
each end, but otherwise two sets of stiffness values must be determined for each 
member, one set for each end. 

The application of the numerical procedure can be outlined in the following 
steps: 


1. The fixed-end moments, about the principal axes, are computed for each 
end of all loaded members. 

2. These moments are resolved to fixed-end moments about the rectangular 
axes of the selected reference system. In general, this will result in unbalanced 
moments about the reference axes at each joint. 

3. At each joint, the unbalanced moment about each axis is distributed, in 
turn, to all members framing into the joint in proportion to their relative stiff- 
nesses about each of the axes. The distribution at each joint is algebraic, and 
signs are determined automatically. 

4. Each distributed moment is multiplied (with due regard to signs), by 
the appropriate carry-over factors in turn, and the resulting moments are car- 
ried over to the proper axes at both the near and far ends. 

5. Steps 3 and 4 are repeated until a desired convergence is reached. 

6. All values following Step 2 are totalled algebraically to obtain the end 
moments about the reference axes. 

7. The values obtained in Step 6 are resolved into bending moments about 
the principal axes and into torsional moments. 

The girder shown in Fig. 7 is fixed at ends A and C, and is continuous over 
support B. Members AB and BC, both of which lie in the xz-plane, are identi- 
cal in length and cross section, with relative bending stiffness about axis 1-1 
equal to ten, and relative torsional stiffness equal to one. Since the load is ap- 
plied normal to the rz-plane, and the angle y is equal to zero, there will be no 
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moments about the y-axis. Asa result of this, it is not necessary to determine 
stiffnesses or carry-over factors associated with that axis. The required con- 
stants are readily determined from the simplified expressions (W = 0;¢ = 90°) 
in Fig.6. Ends A and C are fixed, and there will be no carry-over from A tq B 


or from C to B. 


K,=10 
K,=1 2 
SECTION E-E 


Fie. 7.—Continvovus Beam 


TABLE 1.—ANALysis oF THE ContTINUOUS BEAM oF Fia. 7. 


Moments THE 2-AxIs, Moments THE 2-Axis, 
1n Pounp-Feet 1x Pounp-Feetr 


a Joint B Joint C || Joint A 


Member 
BA 


+0.7 +0.235 
Cn = +0.154) C’sr= 
= +0.799| —1.0 
C = +1.199 


F.E.M. —31,250 +31,250 
Dist. —23,906 —7,344 365% —30,472 
COM (2) —3,682 
COM (z) —7,924 —11,897 
Dist +9,101 +2,796 +12,526 +16,137 
COM (z) +1,402 +10,912 
COM (2) +4,196 +6,301 
Dist. —4,820 —1,481 —4,769 —6,143 
COM (z) —742 8! —5,779 
COM (2) —1,598 —2,399 
Dist. +1,835 +564 +2,525 +3,254 
+283 +2,200 
+846 +1,270 
—972 —298 —961 —1,239 
—150 —1,165 
COM (2) —322 —483 
Dist. +369 +114 +509 +656 
COM (z) +57 4298 +442 
COM (2) +171 +256 
Dist. —196 —60 —193 
COM (zr) —30 ~ 235 
COM (2) —64 —97 
Dist. +74 +103 
COM (z) +11 +89 
COM (2) +34 +52 | 
Dist. —40 —12 —39 
COM (z) —6 —48 
COM —13 —20 
Dist. +15 +5 +21 +27 


Total —38,781 +5,693 | —5,603 —75,258|  +17,947 | —17,947| —8,988 


* Abbreviations: Dist. = Distribution; C.O. = Carry-Over; F.E.M. = Fixed-End Moment; COM (z) 
= Carry-Over Moment from the z-Axis; COM (z) = Carry-Over Moment from the z-Axis. 
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K,=7.75 
y K,=3.25 
q 
Joint B Joint C 
Member Member Member| Member|| Member | Member | Member 
AB BA BC CB AB BC CB 
| Factor: 
Dist. +0.437 +0.563 
C.0. C'ss = +0.335 +0.5 
a C as = +0.503 
~15,236 
+8,068 
q 
‘ —3,072 
+1,627 
—620 
+328 
—25 
: 
| 
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It can be seen by inspection that for member BC, K, is equal to k,:, and 
K, isequal to k,. Likewise, the carry-over factors from B to C are C,, (equal to 
—1.0) and C’,, (equal to +0.5). For convenience in distributing unbalanced 
moments, the ratio of the stiffness of each member to the sums of the stiffnesses 
of the members at each joint, for a particular axis, is entered in Table 1 as a 


distribution factor. For example, the distribution factors for moment about 
3.25 1 


the z-axis at joint B are ——~—— = 0.765 for member BA, and -==—— 


25 


3.2 3.25 + 1 
= ().235 for member BC. 

As previously mentioned, a prime symbol used in connection with a carry- 
over factor indicates carry-over to the far end, whereas the absence of such a 
mark indicates carry-over to an axis at the end about which the distribution was 
made. That is, a moment distributed to an z-axis is multiplied by C,, to ob- 
tain the moment induced about the z-axis at the same end of the member, but 
the distributed moment is multiplied by C’,, to obtain the induced moment 
about the z-axis at the opposite end of the member. 

To determine the fixed-end moments about the x-axis and the z-axis, at both 
joint A and joint B, one must first obtain the fixed-end moments about axis 
1-1. These values are equal to (}) (20,000) 25, or 62,500 lb-ft, at each end. 
The fixed-end moments about the reference axes are 62,500 sin 30° = 31,250 
lb-ft, about the z-axis and 62,500 cos 30° = 54,125 lb-ft, about the z-axis. 
These values, with appropriate signs, are entered in Table 1. 

The tabulated fixed-end moments create a statical unbalance at joint B, 
and the unbalance about each axis is distributed to the two members at B in 
proportion to their relative stiffnesses about the corresponding axes. For ex- 
ample, the unbalanced moment of +31,250 lb-ft at joint B is distributed so that 
23,906 lb-ft is resisted by member BA and 7,344 lb-ft is resisted by member BC. 
The distributed moments are entered in Table 1 with signs opposite to that of 
the unbalanced moment, and the sum of +31,250 lb-ft and —23,906 lb-ft is in 
balance with —7,344 lb-ft. Since joint A is fixed, there can be no unbalance at 
that joint; therefore, no moment is distributed to member AB at that end. 

The carry-over moments (COM) are obtained by multiplying the distributed 
moments by the appropriate carry-over factors. Carry-over moments induced 
by a moment distributed to the z-axis are placed opposite the designation COM 
(x), and carry-over moments induced by a moment distributed to the z-axis 
are placed opposite the designation COM (z). To illustrate, the carry-over 
moments induced by the distributed moment of — 23,906 lb-ft are the following: 
At the far end, joint A, the carry-over moments are COM(x) = — 23,906 
(4-0.154) = — 3,682 lb-ft about the z-axis, and COM (x) = — 23,906 (+0.799) 
= — 19,101 lb-ft about the z-axis; at the near ead, B, COM (r) = — 23,906 
(+1.199) = — 28,663 lb-ft about the z-axis. The distributed moment of 
—7,344 lb-ft induces a moment at the far end, C, whose value is COM (z) 
= — 7,344 (—1.0) = +7,344 lb-ft about the z-axis. 

The carry-over moments resulting from the moments distributed to the 
z-axis at B are obtained in a similar manner, and are recorded opposite the 
designation COM (2). 
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At each joint, the new unbalance resulting from moments carried over is 
distributed in proportion to the relative stiffnesses, and the resulting values are 
entered in the tabulation opposite the designation ‘‘Dist.”” No distribution is 
made at the fixed ends. The procedure of successively carrying over and dis- 
tributing is continued to a desired degree of convergence, and the totals are ob- 
tained. In this example the process has been terminated after a distribution, 
It is a matter of academic interest that the moments balance at joint B, but the 
moments about the z-axis at the ends of member BC are not exactly equal 
(5,693 lb-ft as against 5,698 lb-ft), as torsional moments should be. If the last 
distribution is omitted, the latter moments will balance exactly, but the mo- 
ments at joint B will not. 

Diagrams of bending moments and torsional moments are shown in Fig, 
8(a). The bending moment diagram for AB was drawn by superimposing 
moments resulting from continuity on the moment diagram for beam AB con- 
sidered as simply supported. Since member BC is parallel to the z-axis, the 
bending moments and torsional moments for that member are the moments 


_ Note 


sin®* Os \& All moments are in foot-pounds 
4 
GJ 


4 @ 
\% 
% 
ye 


Torsional! Moments 


(a) Moment Diagrams Plotted Along the Structure 


Fx IF, 


#(125,000)(25)(0.5}(0.5) © *781,250 = <1 1 
388)(25)(0.5)(0.5) -114,925  +(18,388)(25)(0,5)(0.866) = 199,050 
947)(25)(0.5) = 224,337 
~(8,988)(25)(0.5) = 112-350 


EM, (about c) 


-(8,968)(25)(0.5)(8.33) 
(17, 947)(25)(0.5)(16.67) = 3,739,706 
+(18, 388)(25)(0.5)(29.98) = 

(125,000) (25)(0.5)( 3015) = =53,359, 375 
(4, 045)(25)(2.5)(12.5 6 


(b) Numerical Check 


Fig. 8.—Moments and Numerical Check for the Continuous Beam of Fig. 7 
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about the z-axis and the z-axis, respectively, as determined in Table 1. To ob- 
tain the bending moments and torsional moments for member AB, the moments 
obtained in Table 1 were resolved into components about the bending axis and 
the axis of twist. 

A solution of a statically indeterminate structure is correct providing the 
requirements of continuity are satisfied. To satisfy these requirements, the 
angle changes—as determined by means of the moments and the defined proper- 
ties of the material—must “balance” around a closed circuit. By ‘‘balance’”’ is 
meant that the distribution of angle changes must be such that the sum of the 
angle changes about any axis equals zero, and the sum of the products of angle 
changes times distances about any axis at a support—that is, the displacement 


SECTION E-E 


Fie. 9.—Continvous Frame 


TABLE 2.—AnaA OF THE 


Moments Apovut THE z-Axis, IN 


Descriprion Joint A Joint B Joint C Joint D 


Member Member Member Member Member Member 
AB BA BC BD CB DB 


Factors: 
Dist. 0.433 0.134 0.433 
C's = +0.154 C's =-10 C's = +0 500 
Cg = +0.799 = +0.067 
C a = +1.199 C ez = +0.133 


CO. 


F.E.M. —31,250 +31,250 

Dist. — 13,531 —13,531 
COM (2) — 2,084 

COM (2) —6,527 —9,801 1,089 
Dist. +4,715 +4,715 
COM(z) +726 

COM (2) +2,174 +3,264 +363 
Dist —1,570 —1,570 
COM (2) —241 

COM (2) —757 ~-1,137 —126 
Dist. +547 +547 
COM (2) +84 

COM (2) +252 +379 +42 
Dist. —182 —182 
COM (2) —28 

COM (2) 88 —132 —15 
Dist. +64 +64 
COM +10 

COM (2) +29 +44 +5 
Dist —21 


Totals —37,700 | +13,889 | —3,091 ~10,798 +3,084 | —5,392 
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Bending moments plotted 
on tension side 


All moments in foot - pounds, 
but not to scale. 


Torsional Moments Moments About Asis 2-2 


Fria. 10.—Moment Diacrams For THE ContINvovus Frame or Fia. 9 


Continuous FRAME OF Fia. 9 


Moments Asout THE 2-Axts, IN 


Joint A Joint B ‘ Joint C Joint D 


Member Member Member Member Member Member 
AB BA BC BD CB DB 


0.360 0.465 0.175 
Cu= +0452 Cn = +05 Cu= +0.500 
= +0.335 = +0.058 
C = +0.503 +0115 

+54,125 
— 19,485 5 —9,472 
—16,224 — 1,800 


+6,489 ‘ +3,154 
+3,767 +5,653 +627 


+2,933 
—2,261 — 1,099 
—1,254 — 1,882 —209 
— 1,022 
+753 ‘ +366 
+437 +656 +73 


+340 
—262 —128 
—145 —218 


—118 
+87 +42 
+51 +72 
+40 


+27,474 


g 
Nn 
\ 
to 
Note 
=| 
. 
A D & 
Line 
: 
f 
3 —54,125 
5 —10,811 —907 
—8,807 ~ 12,584 —4,736 
2 +4,190 41,577 
11 —105 
12 —1,460 —550 
1 
15 +486 +183 
16 
17 -12 
18 —170 —64 
19 
20 +4 
21 +56 +21 
22 —38 —14 
23 —68,714 — 18,999 ~8,475 —9,482 —4,237 
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at a support—equals zero. Stated algebraically, 


= 5, = + Dydd, = 
5, = + (7) 
5, = + Dydd, = 


In Eqs. 7, d6,, d6,, and d6, are angle changes about the z-axis, the y-axis, and the 
z-axis, respectively ; and 6,, 6,, and 4, are displacements of one end of a circuit 
relative to the other end along the z-axis, the y-axis, and the z-axis, respectively. 

If angle changes are treated as forces along their axes of rotation, a displace- 
ment can be computed as a moment of these “forces” about the axis of dis- 
placement. The requirements of geometry stated in Eqs. 5 can be written as 
analogous requirements of statics, as follows: 


For rotations: For displacements: 


=F, = = 


=F, = 

In Fig. 8(b) a check is made of the solution through application of the ap- 
propriate relationships from Eqs. 8. Angle changes were first obtained about 
the principal axis of bending and the axis of twist using the moments in Fig. 
8(a) and the appropriate values of bending and torsional stiffness per unit 
length of a member. Such angle changes are obtained by dividing bending 
moments by FE J and torsional moments by GJ. It is only necessary to use 
relative values of EIT and GJ. The angle changes about the z-axis and the 
z-axis (“‘forces” in the z and z directions, respectively) were obtained by use 
of the appropriate functions of angle a. In Fig. 8(6), all values of angles are 
based on relative values of E J and GJ equal to 2.5 |b-sq ft and 1.0 lb-sq ft, 
respectively. The results of the check in 8(b) show that the requirements of 
geometry are satisfied. 

Except for the addition of column BD, the structure shown in Fig. 9 and 
analyzed in Table 2 is identical to that in Fig. 7. The column is so oriented 


that the principal axis 2-2 of the section is paralle? with the longitudinal axis of 
member AB. From the expressions in Fig. 6 (¢ = 0), the stiffnesses and carry- 
over factors for BD were determined. The numerical procedure follows in a 
manner similar to that in Table 1. Moment diagrams are shown in Fig. 10. 
A check of the angle changes along the circuit ABC and along the circuit DBC 
showed an almost perfect balance. 

In the frame in Fig. 11 (a), member AB does not lie in any of the three orthog- 
onal planes. The projected angle a of this member is equal to 30°, the angle 
¢ is equal to 105°, and the angle y is equal to zero. Column BD is so oriented 
that the principal axes of its cross section are parallel to the chosen reference 
axes, thus making angle a for this member equal to zero. Girders AB and BC 
are identical in cross section. Values of stiffnesses and carry-over factors for 
member AB were determined through use of the general expressions in Fig. 4. 
Because the principal axes of the cross sections of members BC and BD are 
parallel to the reference axes, their stiffnesses about the latter axes are identical 
to those about the corresponding principal axes. The structure is analyzed in 


Table 3. 
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Note: Support B is Assumed 
to Settle One Inch E = 432000000 Lb per Sq Ft 


1,=0.1447 Ft* 


K,=10 
K,=1 


SECTION E-E 
Fig. 12.—Errecrs or TRANSLATION 


TABLE 4.—ANALYSIS OF THE EFFECTS OF TRANSLATION OF THE 
Continuous Beam oF Fia. 12 


Moments THE z-Axis, Moments ApouTt THE 2-AxI4, 
In Pounn-Feet Pounp-Feer 


D 
von | Joint A Joint B Joint C || Joint A Joint B Joint C 


0.235 
= +0.154| = +0452) = 
= +0.799| —1.0 = +0.335| +0.5 
C az = +1.199 C = +0.§03 


—43,300 


+19,125 
—5,875 || +15,281 


COM (a) —1,473 —1,323 


ist. +1,127 —12,910 
COM (z) —346 +900 
COM (2) —5,041 —4,529 
Dist. +3,856 —761 
COM (z) —1,185 || +3,081 
COM (2) —297 —267 
Dist. +227 —2,642 
COM (z) +181 
COM (z) —1,032 —927 
Dist. +789 —153 
COM (z) +630 
COM (2) —60 —54 
Dist. +46 —533 
COM (z) +37 
COM (2) —208 — 187 

Dist. +159 —31 
COM (z) 4-127 
COM (2) 8 —12 -1i1 
Dist. +9 —108 


Totals —7,785 —7,782 || —30,361 — 29,093 +29,090 


ip 
| 
| 
B cl 2 
1 1 10 
30°~ K,=7.75 2 
Member Member Member| Member|; Member Member Member | Member a = 
AB BA BC CB AB BA BC CB 
Factor: 
Dist. 
| —25,000| —25,000 || —43,300 +50,000 | +50,000 
Dist. —2,928 —3,772 
COM (z) +2,945 +22,931 
f — 1,886 
—380 
—1,321 
— 266 
—16 
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Diagrams of bending moments and torsional moments for the structure are 
shown in Fig. 11. Again a check of the angle changes along the circuits ABC 
and DBC showed a very good balance. Of course, in this case it was necessary 
to check the balance for angle changes (‘‘forces’’) associated with the z, z, and 
y directions. 

The analysis of a structure in which joint translation occurs is illustrated in 
Fig. 12. The structure considered is that of Fig. 7, and it is assumed that the 
support B settles 1 in. Fixed-end moments, about axis 1-1, are equal to 
6 EI = 6 x 432,000,000 x 0.1447 x = 50,000 Ib-ft. Fixed-end 
moments about the reference axes were next obtained and entered in Table 4, 
and the procedure was applied as before. 


SuMMARY 


A procedure has been presented for the structural analysis of girders and 
frames that are continuoys in space. As herein presented, it is applicable to 
structures composed of straight members; the members may be of either con- 
stant or variable cross section. The method involves successive corrections of 
an assumed solution, and consequently lends itself to approximate or precise 


analyses. 
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